On the transmission of light through a single rectangular hole 
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In this Letter we show that a single rectangular hole exhibits transmission resonances that ap- 
pear near the cutoff wavelength of the hole waveguide. For light polarized with the electric field 
pointing along the short axis, it is shown that the normalized-to-area transmittance at resonance is 
proportional to the ratio between the long and short sides, and to the dielectric constant inside the 
hole. Importantly, this resonant transmission process is accompanied by a huge enhancement of the 
electric field at both entrance and exit interfaces of the hole. These findings open the possibility of 
using rectangular holes for spectroscopic purposes or for exploring non-linear effects. 
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Since the pioneering work of Ebbesen et al. Q report- 
ing extraordinary optical transmission (EOT) through 
two-dimensional (2D) hole arrays perforated in optically 
thick silver films, the study of the transmission properties 
of subwavelength apertures (holes or slits) has become a 
very active area of research in electromagnetism. In the 
last few years, EOT has been found in other frequency 
regimes, as THz 0,0 and microwave 0. Nowadays there 
is a wide consensus that EOT in hole arrays is linked to 
the excitation of the surface electromagnetic (EM) modes 
UEl that decorate the surfaces of structured metals. 

Very recently, different experiments have focused on 
the influence of hole shape on the transmission properties 
of both 2D hole arrays 0, H, E| an d single subwavelength 
holes [lOj in the optical regime. All these studies show 
strong polarization dependencies in rectangular holes and 
also that rectangular holes exhibit higher transmittance 
than square or circular holes with the same area. Inter- 
estingly, it was also found that single rectangular holes 
can support transmission resonances, even in the sub- 
wavelength regime. 

In this Letter we present the, up to our knowledge, 
first theoretical study on the dependence on hole shape 
of the transmittance through a single hole. We find that, 
as a difference with circular holes [llj , single rectangular 
holes can exhibit strong transmission resonances in all 
frequency ranges. One of these resonances appears close 
to cutoff, with a peak transmittance controlled simply by 
the ratio between the long and short sides of the rectan- 
gle. Additionally, we show that the presence of a material 
filling the hole greatly boosts the transmittance. Asso- 
ciated to these transmission resonances, there is a very 
strong increase of the electric field at the hole. 

Figure I shows schematically the system under study. 
A rectangular hole of sides a x and a y perforated on a 
metallic film of thickness h. The system is illuminated 
by a p-polarized plane wave with wavelength A, the in- 
planc component of the electric field pointing along the x- 
direction. The metal is treated within the perfect conduc- 
tor approximation (PCA). This approximation is known 
to fail in two instances: first, when the dimensions of 




FIG. 1: Diagram of a single rectangular hole of sides a x and 
a y perforated on a metal film of thickness h. The structure 
is illuminated by a p-polarized plane wave with the magnetic 
field pointing along the y-direction. 



the structure are of the order or smaller than the skin 
depth [l^. Even in this case, the range of validity of 
the PCA can be greatly extended by simply consider- 
ing effective hole dimensions enlarged by the (metal and 
wavelength dependent) skin depth 0| . Secondly, in the 
optical regime there are effects related to the existence of 
absorption and surface plasmons, which are not captured 
by the PCA. However, the PCA serves both as the start- 
ing point for more elaborated approximations, (as the 
one considering surface impedance boundary conditions) 
and for clarifying which effects come from geometry and 
which have a dielectric origin. With these caveats in 
mind, especially important in the optical regime, the re- 
sults presented in this letter apply to different frequency 
regimes, by simply re-scaling all lengths by the same fac- 
tor. 

Let us briefly describe the formalism used for calcu- 
lating the transmittance through the structure (a de- 
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tailed account of this method that was developed to treat 
an arbitrary number of indentations can be found in 
[f|). In this method, the EM fields in both reflection 
(I) and transmission (III) vacuum regions are expressed 
in terms of EM eigenmodes \kcr >, characterized by the 
in-plane component of the wavevector k, and the polar- 
ization a. Inside the hole, the EM field is expanded in 
terms of all EM waveguide eigenmodes. After match- 
ing appropriately the EM fields at the two interfaces of 
the system (z = and z = h), the formalism provides 
the full EM field in all space in terms of the projection 
onto waveguide eigenmodes of the electric field at both 
hole entrance and exit interfaces. In all calculations pre- 
sented in this letter (normal incidence illumination and 
a x ,a y < A), we have checked that considering just the 
first TE eigenmode (\TE >) is enough to obtain very 
accurate results for the transmittance so, for simplicity, 
we present our formalism just for this case. Then, in 
terms of the modal amplitudes E and E' , the electric 
field bivector E = (E x ,E y ) T (T standing for transpo- 
sition) at the hole entrance and exit can be written as 
\E(z = 0) >= E \TE > and \E(z = h) >= -E' \TE >, 
respectively. Here we have used Dirac's notation, with 
< r \TE >= (1, 0) T sin[7r(y/a a + l/2)]/vW, Af = a x a v /2 
being a normalization factor. The equations that E and 
E' must satisfy are: 

{G — Yi)E — G v E' = I , 
-G V E + (G-Z)E' = 0. (1) 

where Iq — 2i < kop \TE > is the external illumination 
term, i.e., the overlap integral between the incident plane 
wave \kop > and mode \TE >. For the case of normal 
incidence, and normalizing the incident EM field such 
that the incoming energy flux over the hole area is unity, a 
simple calculation gives Iq = 4v / 2«/ 7r - ^ an d Gv are two 
magnitudes that only depend on the characteristics of 
TE mode inside the hole: S = Yte/ ta,n(q z ti) and Gy = 
YteI sin(q z h); q z = yjk 2 ~ ( n / a y) 2 i s tne propagation 
constant of the fundamental TE mode, Yte — Qz/k^ its 
admittance and k u = 2ir/\. 

The self-illumination of the hole, via vacuum modes, 
is controlled by G = i/(2ir) 2 I d ^ Y kJ < TE \k<J > | 2 
(with Yg being the admittance of the EM vacuum mode 

\ka > 15] ). For the case of rectangular holes, 

x [ Sinc( h^y±JL) + sinc ( V|Z^)]2 (2 ) 

where k 2 = k 2 + k 2 . Notice that, in our formulation, 
Rc(G) comes from the coupling to evanescent modes in 
vacuum and Im(G) from the radiative modes. 

After obtaining E and E' from (1), the normalized-to- 
area transmittance (T) is calculated as: 



T=^lm[E*E>] (3) 

k p 

Figures 2 and 3 illustrate the dependence of T with 
the two geometrical parameters involved (ratio a y /a x and 
thickness h). In both cases we consider normal incidence. 
As previously said, all results are scalable, and we have 
chosen as unit length the cutoff wavelength of the TE 
waveguide inside the hole, A c = 2a y . 




Wavelength (in units of X c ) 

FIG. 2: Normalized-to-area transmittance (T) versus wave- 
length (in units of the cutoff wavelength, A c ), for a normal 
incident plane wave impinging at a rectangular hole, for dif- 
ferent ratios a y /a x . Metal thickness is h — a y /3. For com- 
parison, in the inset we plot T versus wavelength for a single 
square (black line) and circular (red line) holes. 

Figure 2 renders T for the case h = a y /3, for different 
values of a y /a x . As clearly seen in this figure, a trans- 
mission peak develops at approximately A c , with increas- 
ing maximum transmittance and decreasing linewidth as 
a y /a x increases. In the case of square or circular holes 
there is also a resonance close to cutoff, but a very faint 
one (see inset of Fig. 2). In these last two cases, bellow 
cutoff the normalized-to-area maximum transmittance is 
of order of 1, i.e., approximately the amount of light that 
is directly impinging at the hole. In all cases, above cut- 
off T decreases strongly with A, both due to the fact that 
the fields inside the hole are evanescent and that, in the 
extreme subwavelength regime, an incident wave couples 
very poorly to the hole [ljj ■ 

The dependence of T(A) with metal thickness is shown 
in Fig. 3, for a rectangular hole with a y /a x = 10. Apart 
from the previously discussed transmission peak located 
at A ~ A c , which spectral position essentially remains 
independent of h, a series of transmission resonances 
emerge as the depth of the hole is increased. As we 
will show later on, these are Fabry-Perot resonances in 
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FIG. 3: T for a normal incident p-polarized plane wave versus 
wavelength (in units of A c = 2a y ) and thickness h (ranging 
from two 4a y ) for rectangular holes with aspect ratio a y = 
10a x . 



rectangular holes, similar to the ones recently found in 
sub- wavelength ID slits [HI El El 

The resonant characteristics of the transmittance 
through rectangular holes, and their dependence on geo- 
metrical parameters can be worked out analytically from 
set of equations (1). For the case we are analyzing (a sym- 
metric structure with respect to the plane z = h/2), max- 
imum transmission appears when the field intensities at 
the entrance and exit sides of the aperture are equal, i.e. 
\E\ = \E'\. From (1), this occurs when \G - S| = \G V \, 
condition that after some straightforward algebra implies 



2Re(G) = J-L; ^tan(g z /i) 



Yn 



TE 



(4) 



There are several wavelengths at which this transcen- 
dent equation is satisfied, providing the exact spectral 
location of the transmission peaks (A res ) found in Figs. 
2 and 3. Ignoring the shift in the spectral dependence 
due to the EM coupling to vacuum modes (this is, setting 
G — > 0, which is the appropriate limit for a x ,a y << A), 
Eq. (4) transforms into Yte t&n(q z h) = 0, which is the 
usual Fabry-Perot condition for the existence of a stand- 
ing wave inside the hole. Note that this last equation al- 
lows the solution q z — 0, so a transmission peak located 
at around the cutoff wavelength is expected, irrespective 
of the geometrical parameters a y /a x (see Fig. 2) and h 
(see Fig. 3). In general, Eq. (4) predicts the shifts in the 
Fabry-Perot resonant wavelengths due to the coupling to 
both radiative and evanescent vacuum modes. 

Using the resonance condition [Eq.(4)], from Eqs.(l) 
and (3) and after some straightforward algebra, we obtain 
that the normal-incidence T at resonance, T res is given 
by: 



T — 



4Im(G) 



(5) 



A very good simple approximation to T res can be ob- 
tained realizing that, in the extreme subwavelength limit 
(a x ,a y « A), Eq. (2) gives Im(G) « 32^^/(3^ A 2 ). 
We have checked that this expression holds even for 
flj a: A/2 therefore, for A > 2a y we find 



T 



3 A 2 
4-7T a x dy 



(6) 



Recalling that T is the normalized-to-area transmis- 
sion, this expression implies that the total amount of 
light emerging from a rectangular hole is, at least for 
the resonance appearing close to cutoff, independent of 
the length of the short side! Although derived for rect- 
angular holes, Eq. (6) seems to be more general as the 
same expression was found for circular holes 201, with 
the term a x a y replaced by the area of the circular hole. 
The important point in rectangular holes is that, for the 
polarization chosen, the transmittance peak appearing at 
cutoff depends only on the long side (A res = 2et y ), result- 
ing in a transmittance T res « (3/Tr)a y /a x close to cutoff. 
This is the main result of this Letter, as it predicts a 
huge transmission enhancement in a single rectangular 
hole with large aspect ratio. 




Wavelength (in units of 2a ) 



FIG. 4: T for a normal incident plane wave versus wavelength 
for a rectangular hole with a y /a x — 10 and different values of 
e inside the hole. Metal thickness is ft = a y /3. Dashed and 
dotted lines show the behavior of Eqs.(5) and (6), respectively. 
Orange curve shows T versus A in the limit h — > for a 
rectangle with aspect ratio a y /a x = 10. 

Additionally, even for a fixed aspect ratio a y /a x , Eq. 
(6) gives us a clue for further enhancing the transmission: 
filling the hole with a material with dielectric constant 
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e > 1, as this increases the cutoff waveiength. In fact, 
in the definition of quantities appearing in Eqs. (1) and 
(3), the only place in which e enters is in the propagation 
constant associated to mode \TE > which now reads: 
q z = Weft* — (ir/ay) 2 . Therefore, the spectral position 
of resonances depend on q z (and therefore on e) but the 
transmittance at resonance is still given by Eq. (5). As a 
result, filling the hole with a dielectric would redshift the 
Fabry-Perot transmission peak appearing close to cutoff 
to as 2y r ea y (for thick enough metal films, see below) 
and, more importantly, increase its transmittance. This 
is illustrated in Fig. 4, which renders the transmission 
spectra for rectangular holes of aspect ratio a y /a x = 10, 
in a metallic film of thickness h — a y /3, for several values 
of e. Note that this way of increasing the transmission 
through the hole by filling it with material with e > 1 can 
be also operative for the case of circular [2l| or square 
holes. Remarkably, in rectangular holes this mechanism 
acts almost independently of the enhancement due to the 
aspect ratio, so T res is proportional to both a v /a x and e. 

It is worth analyzing the limit h — > of the transmis- 
sivity of single rectangular holes, the analogue of Bethe's 
study for circular holes [16j. By taking this limit in the 
set of equations (1) and Eq.(3), we can obtain an expres- 
sion for T valid for all wavelengths: 



These effects (huge enhancements of both transmission 
and field amplitude) should be readily observable in the 
microwave or THz frequency ranges. In these regimes, 
holes with very large aspect ratio can be manufactured 
and, furthermore, dielectric materials presenting large 
positive dielectric functions are also available. 

We would like to end with some comments about the 
transferability of our results to the optical regime. Al- 
though our results agree with the experimental finding 
that transmittance resonances increase and become bet- 
ter defined with increased aspect ratio, they have only 
a semi-quantitative value in this regime, especially for 
short sides not much larger than twice the skin depth 
(i.e. « 50nm). Additionally, the influence of localized 
surface plasmon modes on the transmittance is not cap- 
ture by our model. The role played by these modes on the 
transmittance is a point that deserves further theoretical 
investigation. 
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4\G\ 



r Im(G), 



(7) 



which, as expected, is independent of e. The behavior of 
this magnitude as a function of wavelength is rendered in 
Fig. 4 (orange curve) for the case a y = 10a x . As metal 
thickness is decreased, the resonance moves from a loca- 
tion close to 2y/ea y to a much shorter wavelength, given 
by the resonant condition Re(G) = (notice that Eq.(7) 
at resonance gives the same expression as Eq.(5) for this 
case). Interestingly, in the extreme subwavelength limit, 
as Re(G) oc A and Im(G) cx 1/A 2 , T decreases with A as 
1/A 4 , as in the case of circular holes (lif . 

It is also interesting to analyze the enhancement of 
the EM-fields associated to this resonant phenomenon. 
Naively, one would expect that the intensity of the E- 
field at the entrance and exit sides of the hole (\E\ 2 and 
\E'\ 2 ) should be proportional to the transmittance. How- 
ever, the direct evaluation of \E\ 2 = \E'\ 2 at the resonant 
condition given by Eq. (4) yields: 



\E\: 



\M 



(8) 



4[Im(G)] 2 ' 

leading to an enhancement of the intensity of the E-field 
(with respect to the incident one) that scales with A res 
as X^ es /(a y a x ) 2 , square of the enhancement in the trans- 
mittance (see Eq.(6)). This implies that in the process 
of resonant transmission, light is highly concentrated on 
the entrance and exit sides of the hole but only a small 
fraction of this light is finally transmitted. 
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